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' We compare the apparent axial ratio distributions of Brightest Cluster Galax- 

ies (BCGs) and normal ellipticals (Es) in our sample of 75 galaxy clusters from the 
WINGS survey. Most BCGs in our clusters (69%) are classified as cD galaxies. The 
sample of cDs has been completed by 14 additional cDs (non-BCGs) we found in our 
clusters. We deproject the apparent axial ratio distributions of Es, BCGs and cDs using 
a bi-variate version of the rectification Lucy's algorithm, whose results are supported 
by an independent Monte-Carlo technique. Finally, we compare the intrinsic shape 
distribution of BCGs to the corresponding shape distribution of the central part of 
cluster-sized dark-matter halos extracted from the GIF2 ACDM iV-body simulations. 

We find that: (i) Es have triaxial shape, the triaxiality sharing almost evenly 
the intrinsic axial ratios parameter space, with a weak preference for prolateness; (ii) 
the BCGs have triaxial shape as well. However, their tendence towards prolateness 
is much stronger than in the case of Es. Such a strong prolateness appears entirely 
due to the sizeable (dominant) component of cDs inside the WINGS sample of BCGs. 
In fact, while the 'normal' (non-cD) BCGs do not differ from Es, as far as the shape 
distribution is concerned, the axial ratio distribution of BCG_cD galaxies is found to 
support quite prolate shapes; (iii) our result turns out to b e strongly at variance with 
the only similar, previous analysis by iRvden et al.l (|l993l RLP93), where BCGs and 
Es were found to share the same axial ratio distribution; (iv) our data suggest that the 
above discrepancy is mainly caused by the different criteria that RLP93 and ourselves 
use to select the cluster samples, coupled with a preference of cDs to reside in powerful 
X-ray emitting clusters; (v) the GIF2 TV-body results suggest that the prolateness of 
the BCGs (in particular the cDs) could reflect the shape of the associated dark matter 
halos. 

Key words: galaxies: clusters - galaxies: general - galaxies: elliptical and lenticular 
- galaxiesxD 



1 INTRODUCTION 
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Rich galaxy clusters show in their central part a remark- 
able concentration of galaxies, surrounded by progressively 
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less dense regions. The brightest cluster galaxies (BCGs) 
are usually elliptical-like galaxies, often much brighter than 
the rest of the global popula tion of Es ([Sandage and Hardvl 
1 19731 ; ISchombertll 19861 . 119921 ) . They are nor mally found close 
to th e peaks of galaxy nu mber density dBeers and Gellerl 
Il983h and X-ray emission ll Jones and Formanl 11984 1 of the 
cluster. The significant alignment between the elongations 
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with the correlations between BC Gs luminosity and cluste r 
properties (i.e. X-rav temperature; lEdge and Stewartlll99ll ). 
suggests that the formatio n history of BCGs is closely linked 
to that of the host clusters l|Kormendv and Diorgovskill 19891 ; 
ICoziol et alj|2009h . 

Among the different formation scenarios which have 
been proposed in the literature for BCGs, we mention: 
(i) merging of compact galaxy gr oups before or during 
clust er assembly a n d virialization |Ostriker and Tremaind 
1 19751 ; IWhitd 1 19761 ; iMerrittl 1 19851 . see also TV-body sim- 
ulations by Dubinski 1998); (ii) filament driven, regu- 
lar a ccretion of small, gas-rich proto-galactic units |Westl 
1 19941); (iii) tidal stripping f rom other clus t er ga laxies 
^Gallagher and Ostrikerl 1 19721 ; iRichstoiid" Il975l . Il976h ; (iv) 



continuous accretion of star for ming gas f rom the intra - 
cluster medium (cooling flows; ISilkl Il976l ; iFabianl 1 1994) ; 
(v) late assembling (z ~0.5) of smaller, gas-poor (red) 
galaxies in hierarchic al scenario l|De Lucia and BlaizoifeOQTl ; 
iBernardi etafl l2007h . Although each one of the proposed 
mechanisms turns out to present some drawbacks when com- 
pared with the observations, all of them are likely to con- 
tribute, with different strength and timings, to the whole 
process of formation of BCGs. 

From the observational point of view, the BCGs are spe- 
cial in many aspects: besides the unusually high luminosity, 
the BCGs exhibit huge sizes, sometimes further enlarged by 
diffuse, very extended halos. In this case they are c alled cD 
galaxies and their size s may even reach ~300 kpc (|Oemlerl 
Il976t ISchomberj|l988h . The BCGs are also extremely mas- 
sive objects (up to ~ 1O 13 M0), with peculiar kinematics, in 
that they have lower velocity dispersions and larger radii 
than pr edicted by the Faber- Jackson and Kormendy re- 
lations dThuan and Romanishin] Il98ll; [ Hoessel et al.l 1 19871; 
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20071 ). consistent with the presence of a larger fraction 
of dark matter (|von der Linden et al.l 120071 ) and/or with 
significant growth of BCGs via dissipationless mergers 
|Desroches et alj |2007m. Finally, the BCGs often di splay 



multiple nuclei (ISchneider et al 1 19831 ; lLaine et al.ll2003f ) and 
have been frequently id entified as powerful radio sources 
|Giacintucci et al.ll2007l ). 

The fact that BCGs are peculiar in so many aspects 
with respect to the population of normal Es, might suggest 
that their shape too differs from that of Es. Also, the spe- 
cial role played by BCGs in the formation history of clusters 
could result in a rather peculiar shape. Normal ellipticals 
are g enerally believed to cover the whole range of triaxi- 
ality (|Fasano and Vid Il99ll ; iRvdenl 1 19921 : iBak and Statlerl 
2000), while both evolutionar y scenarios supported by sim- 
ple dynamical considerations (West 1994) and TV-body simu- 
lations of self-consistent models of gala xy clusters including 
dark matter component (|Dubinskill 998). both predict BCGs 



with significantly prolate shapes. Remarkably enough, also 
the TV-body simulations of BCGs-scale, dark-matter ha- 
los in ACD M cosmology lead to similar conclusions about 
their shape (IWarren et al.l 19921; iBailin and Steinmetzll2005l ; 



lAllgood et al2006l : Bett et al 1l2007f l. with the additional, in- 
teresting hint that the dark matter halos become more pro- 
late towards their inner part and at increasing the halo mass 
l|Cole and Lacevlll996l ; iGottlober and Yepedl2007h . In spite 
of these converging indications, the only direct comparison 
available up to now in the literature bet ween the observed 
axial ratio distributions of BCGs and Es |Rvden et al .11 19931 . 
hereafter RLP93, observations in the Kron-Cousins R band) 
led to the conclusion that they are very similar. 

In this paper we perform a new comparison between 
the distributions of the observed axial ratios (hereafter 
g=minor-/major-axis) of BCGs and Es exploiting the large 
databas e of nearby cluster galaxies provided b y the WINGS 
survey |Fasano et al.ll2006l ; IVarela et al.ll2009l ). We pay par- 
ticular attention to the distinction between normal BCGs 
and cDs, thus including in the sample the non-BCGs, cD 
galaxies. We also deproject the apparent axial ratio distri- 
butions to get the corresponding distributions of the intrin- 
sic axial ratios. Finally, using the GIF2 ACDM cosmological 
simulations, we extract 3510 cluster-sized dark- matter halos 
and calculate the central shape at the corresponding BCG 
scale, allowing us to compare their intrinsic axial ratio dis- 
tributions with those of the observed BCGs. 

In Section [2] we describe the galaxy sample and the ax- 
ial ratio data we use in our analysis. In Section [3] the axial 
ratio distributions of Es, BCGs and cDs are presented and 
discussed and the comparison with the results of RLP93 is 
performed. Section [4] outlines the technical aspects of the 
deprojection and presents the distributions of the intrinsic 
axial ratios we obtain applying such deprojection techniques 
to the Es, BCGs and cDs samples. In Section [5] we discuss 
our results and compare the intrinsic shapes of BCGs to 
those obtained from the GIF2 TV-body data. Finally, Sec- 
tion[6]summarize our conclusions. Throughout the paper we 
use the following cosmology: Ho=70 kms _1 Mpc _1 , Qm~0.3 
and Q A =0.7. 



2 GALAXY SAMPLES AND AXIAL RATIO 
DATA 

We have extracted our samples of BCGs and Es from 
the Wide-field N earby Galaxy-cluster Survey (WINGS; 
iFasano et all 120061 ). The 77 clusters of the WINGS sample 
were selected, in the redshift range 0.04-0 . 07, from RO SAT 
catalogs of clusters labeling et alj Il996l . Il998l . |2000| ) and 
turn out to cover a wide range of masses (~ 5 x 10 13 - 
3.2 x 10 15 M ; LogLx ~43.3-44.7 ergs s" 1 ). The opti- 
cal WINGS survey provides B- and V-band imaging of the 
whole sample of clusters. Integrated and aperture photome- 
try hav e been obtained for ~400,0 00 galaxies by using SEx- 
tractor l|Bertin and Arnoutd fl996). The procedures adopted 
to avoid mutual photometric contamination between big 
galaxies with extended halos and sm aller, halo-embedde d 
companions are described in detail in IVarela et all {2009). 
We just recall here that the largest galaxies in each clus- 
ter were carefully modeled with IRAF-ELLIPSE and re- 
moved from the original images in order to allow a reliable 
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masking of the small companions when performing the sur- 
face photometry of the big galaxies themselves. The BCG 
in Abell 3164 has been excluded from the sample since its 
surface photometry turned out to be uncertain due to the 
proximity of the inter-chip region of the CCD mosaic, while 
the BGC in Abell 3562 was not included due to the quality 
of the WINGS imaging for this cluster is not good enough 
to allow a reliable surface photometry and morphology esti- 
mate. Thus, the final sample of BCGs consists of 75 galaxies. 

We have class ified their mo r pholo gy using the WINGS 
V-band imaging |Fasano et al.l 120061 ) and the purposely 
devised, automatic too l MORPHOT (Fasano et al. 2010, 
in preparation; see also 



2009;; IValentinuzzi et al. 



Fasano et al.l 120071 : iPoggianti et al.l 
2009). The logical sequence and the 



basic procedures of MORPHOT are outlined in the Ap- 
pendix, with particular reference to its capability of disen- 
tangling cDs from BCG_Es. 

On the other hand, the quantitative, multi-component 
analysis of the luminosity profiles has been performed, 
on both the V- and B-band WINGS images, using 
GASP2D (Ascaso et al. 2010, in preparation; see also 
iMendez-Abreu et alj [20081 ). In this case, cD galaxies have 
been identified by the simultaneous occurence of two con- 
ditions: (i) the presence, in the outer profiles, of a light 
excess (with respect to the inner Sersic component) signifi- 
cantly larger than the photometric errors; (ii) the positioning 
around a surface brightness level of \iv ~24 of the 'break- 
ing point' where the profile splits up from the inner profile. 
When MORPHOT and GASP2D (in both V- and B-band) 
produce different results or when both tools fail to converge 
(~33% of the BCG sample), we have assigned the morpho- 
logical type (E/cD) relying on the visual inspection of the 
V-band images. The above procedure led us to classify most 
of the BCGs in our clusters (52 galaxies) as cDs (BGC_cD, 
hereafter). Table [1] reports the basic information about the 
BCGs in our sample, as well as some relevant properties of 
the host clusters. 

In the clusters Abell 3395 and Abell 3556 the BCGs are 
quite off-centered with respect to the main concentration of 
galaxies, while the brightest objects belonging to such con- 
centration (in both cases cD galaxies) turn out to be just 
slightly fainter than the corresponding BCGs. This moti- 
vated us to examine the very luminous galaxies in each clus- 
ter, searching for more non-BCG, cD candidates (nBCG_cD, 
hereafter). We found 14 additional cDs (including the two 
cDs previously mentioned in the central part of Abell 3395 
and Abell 3556), all of them being just slightly fainter than 
the corresponding BCGs (see Table [5]). Together with the 52 
BCG_cDs in Table [TJ these objects form a complete sample 
of 66 cD galaxies to be compared with the BCG_E sam- 
ple and with the sample of normal Es extracted from the 
WINGS morphological catalogs (MORPHOT). We decided 
to include in the sample of Es just galaxies with absolute 
V magnitude My <-19.5 and cluster-centric distance less 
than O.6-R200E} The first condition guarantees that the mor- 



1 R200 i s defined as the radius delimiting a sphere with interior 
mean density 200 times the critical density, approximately equal 
to the cluster virial radius. O.6-R200 roughly corresponds to iisoOi 
whose interior mean density is 500 times the critical density. Only 
a few WINGS clusters have pho tometric coverage slightly smaller 
than O.6R200 JCava et al.ll2009h 



phological classification is robust enough, while the second 
condition is necessary to homogenize the photometric cov- 
erage of the clusters. Both conditions make the field con- 
tamination practically negligible. After having removed the 
BCG_Es and the cDs we are left with a final sample of 1024 
Es. 

In the framework of the WINGS survey, the sur- 
face photometry of several hundreds galaxies per clusters 
(those with projected area greater than ~22 arcsec 2 at 
the isophotal level corresponding to 2.5 times the r.m.s. of 
the background) has b een obtained by using GASPHOT 
|Pignatelli et alj 120061 ; iD'Onofrio et all 120091 ). This pur- 
posely devised tool performs simultaneous best fitting of the 
major- and minor-axis growth light curves of galaxies with a 
2D flattened Sersic-law, convolved with the appropriate (lo- 
cal) PSF. In the present paper we use for the shape analysis 
the axial ratios q G of the Sersic model coming out from this 
best-fittig procedure. The reason for this choice is twofold: 
(i) since q G is derived from a best-fitting procedure on the 
whole galaxy body, it is more stable than whatever isophotal 
axial ratio; (ii) since our result is at variance with RLP93 
in that we find more flattened BCGs, and since q G pro- 
vides the lowest flattening among the various possible axial 
ratio estimates (as shown below; see Figure [TJ , our choice 
is a conservative one. Besides q G , GASPHOT provides the 
axial ratio profiles of the galaxy isophotes. In Table [1] and 
[5] we report the isophotal axial ratios of the BCG and cD 
galaxies corresponding to the semi major-axes of 15, 30 and 
60 kpc (q 15 , q 30 and <Z 60 ). It is worth noticing that q G is not 
seeing-affected by definition, since it comes from the best fit 
of PSF-convolved Sersic models. Instead, this is not true for 
the isophotal axial ratios, at least in the innermost galaxy re- 
gions. However, we can assume the seeing influence to vanish 
for D^ oph /FWHM>3, where D^ oph is the isophotal minor- 
axis. In the case of our (BCG_E+cD) sampl e, even in the 
worst seeing conditions (FWHM~2 arcsecs; see lFasano et al.l 
2006) and for the most flattened objects (b/a~0.3), the ax- 
ial ratios of isophotes with semi-major axis >10 arcsecs can 
be considered not seeing-affected. For most distant WINGS 
clusters (z~0.07) this translates into ~13.5 kpc of isophotal 
semi-major axis. We conclude that all the isophotal axial 
ratios reported in Tables [JJ and [5] are not seeing affected. 
We have also collected, for our (BCG-E+cD) sample, the 
axial ratios q s from our photometric catalogs (jVarela et al.l 
I2009I . SExtractor) and, if avail able, the axial ratio s q L from 
the LEDA Hypercat database l|Paturel et alJlioOSl see again 
Tables [JJ and [2} . The last ones are actually literature data 
statistically normalized to the isophote B25, while the first 
ones are luminosity weighted axial ratios similar to those 
defined by RLP93, but referred to the whole galaxy body 
(the axial ratios from RLP93 are instead computed within 
~13 kpc). All these additional estimates of q are useful to 
understand how the axial ratio definition could influence our 
results. 

In Figure [1] the six different estimates of the axial ra- 
tios of our BCG_Es+cDs are compared among each another. 
The BCG_Es, BCG_cDs and nBCG_cDs are represented in 
the figure by white, black and grey (green in the electronic 
version) dots, respectively. For clarity, we will identify the 
panels with the numbers placed in bottom right. First of all, 
from plots 4, 5 and 9 it is clear that in our (BCG_E+cD) 
galaxy sample the flattening of the isophotes usually in- 
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Figure 1. Comparison among different axial ratio measurements for the BCG_Es (white dots), BCG_cDs (black dots) and nBCG_cDs 
(grey dots; green in the electronic version) of the WINGS survey. 



creases towards the outer galaxy regions. From plots 13 and 
15 it is also evident that the axial ratios q G are sistemati- 
cally larger than both q s and q L . From plots 2 and 10 the 
axial ratios q s and q G turn out to be consistent with q 60 
and q 15 , respectively. Finally, the best consistency with q L 
is found, although with a large scatter, in panel 8 with q 30 . 



3 THE AXIAL RATIO DISTRIBUTION OF 
BCG AND NORMAL ELLIPTICALS 

In Figure[2]we compare the axial ratio distributions <f>(q G ) of 
Es and BCGs (panel a), BCG_Es and cDs (panel b), BCG_Es 
and Es (panel c), nBCG_cDs and BCG_cDs (panel d). The 
full line (red in the electronic version of the paper) in the 



top-leftmost panel of the fig ure, reports the a xial r atio dis- 
tribution of Es obtained by iFasano and Viol (|l99ll ). which 
turns out to be fairly in agreement with that of WINGS Es. 
The bottom panels show the corresponding cumulative dis- 
tributions and report, in each case, the probability of the 
"null hypothesis" (i.e. that they are drawn from the same 
parent populations), according to the Kolmogorov-Smirnov 
(KS) statistics. Figure (panel a) shows that the distri- 
bution <f>(q G ) of BCGs looks strongly different from that 
of normal ellipticals, the BCGs being sistematically flatter 
than Es. 

2 Note that our sample of Es includes galaxies with absolute V 
magnitude brighter than -19.5: a rather faint cutoff, indeed. Since 
luminous Es (Mb <-20, i.e. roughly My <-21) are found to be 
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axial ratios from GASPHOT 



Figure 2. Differential (top panels) and cumulative (bottom panels) axial ratio distributions in the WINGS survey: comparisons between 
Es and BCGs (leftmost panels; red and black, respectively), BCG_Es and cD galaxies (panel b), BCG_Es and Es (panel c), nBCG.cDs 
and BCG_cDs (rightmost panels). The bottom panels also report the probability that they are drawn from the same parent populations, 
according to the KS statistics. 



This result is strongly at variance with that obtained 
by RLP93, who conclude that BCGs and Es share the same 
axial ratio distribution. Trying to interpret this quite evi- 
dent discrepancy, we could invoke the different procedures 
adopted to measure the axial ratios. However, we should 
note that, among the six different values of q we have col- 
lected for our BCGs sample, the one we use to get our distri- 
butions (q G ) is by far the one producing the least flattened 
galaxies (see Figure Q]), being also quite consistent with q 15 
(the same outermost isophote used by RLP93). One could 
object that, with respect to our method, the one used by 
RLP93 to measure q tends to overweight the inner galaxy 
regions, which for the BCGs are well known to be rounder 
than the outer ones (see again Figure[T]). However, we should 
also note that our SExtractor estimates, even if computed 
in a fashion similar to RLP93, are actually in fair agreement 
with the isophotal axial ratios at 60 kpes (q 60 ; see plot 2 in 
Figure [l]), which give the flattest distribution. Again, one 
should object that our SExtractor estimates are computed 
on the whole galaxy body, while the RLP93 ones refer to 
the region inside the 15 kpes. In any case, it is worth stress- 
ing that, using the same axial ratio definition (q G ), we get 
very different q distributions for BCGs and normal Es. At 
this point, the discrepancy between our result and that ob- 
tained by RLP93 remains rather puzzling. In Section [5] we 
will present a strong guess about the origin of such discrep- 
ancy. For now we just note that, besides the result concern- 
ing the shape difference between BCGs and Es, the Figure [2] 
is telling us something else: (i) cD galaxies are significantly 

rounder iTremblav and Merritt Il996ft . had we adopted a brighter 
cutoff, the difference between the axial ratio distributions of Es 
and BCGs would have been even larger. 



flatter than BCG_Es (panel b), thus implying that the BCGs 
do not constitute an homogeneous class of objects; (ii) the 
distribution of BCG_Es does not significantly differ from 
that of normal Es (panel c); (iii) the cD galaxies have an 
unique, characteristic shape distribution, regardless of the 
luminosity ranking in the cluster (panel d). Note that this 
result basically justifies why in panel b BCG_Es are com- 
pared with the whole sample of cDs rather than just with 
BCG.cDs. 



4 INTRINSIC SHAPES OF BCG AND 
NORMAL ELLIPTICALS 

4.1 Technicalities 

Following RLP93, hereafter we assume that the intrinsic iso- 
luminosity surfaces of a galaxy are similar coaxial ellipsoids 
with axis lengths in the ratio l:/3:7, with 1JJ ft 7, P and 
7 being the intrinsic axial ratios of the galaxy. Moreover, 
following iFranx et al.l |l99ll ) , we quantify the triaxiality of 
galaxies through the parameter T = (1 — p 2 )/(l — r ) 2 ) (which 
becomes 0/1 for perfectly oblate /prolate bodies) and illus- 
trate the results in the (/3,7) plane (see also Figure 1 in 
iKimm and YHl2007l ). 

To explore the probability density funtion il>(P, 7) of the 
intrinsic axial ratios for a given sample of galaxies (hopefully 
of the same morphological family), we deproject the distri- 
bution of their observed axial rati os 4>(q n ) usin g the iterative 
rectification algorithm devised bv lLucvl (|l974l ). in particular 
its equ ations (13), (14) and (17). In our case the equation 
(17) in lLucvl (l974l ) becomes: 
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r +1 (P,l) = ^J2^^-P[Qn\P,7,<Qn)], (1) 
n — 1 

where N is the size of the galaxy sample and r is the itera- 
tion number. In this equation the r-th estimation of 4>{q) is 
obtained through the equation: 

r (gn) = y J \p r (l3, 1 )P[q n W, 1 ,e(q n )]dpd-f, (2) 

where P[q\/3, 7, e(g)] is the probability that an ellipsoid with 
intrinsic axial ratios /3 and 7 is observed with apparent axial 
ratio q, assuming a random space orientation and the r.m.s. 
of the errors (e) in the q measurements to be a function of 
q itself. 

The iterative machine requires a first guess ^>°(/3, 7) 
of the probability density function we are ex ploring and a 
choice of the error function e(q). According to lFasano et al.l 
i| 19931 . see Figure [T] and equation 1 therein), we assume the 
last one to be a linear function of q. In particular, after some 
statistical comparison among literature data, we decided to 
use the equation: e(q) = 0.05(1 + ?). Concerning the choice of 
t/> (/3, 7), we tried different analytical functions, concluding 
the most flexible and robust guess to be a confined normal 
function, bi-variate in the two quantities 7 and T (triaxi- 
ality). In this way, the first guess of ip requires eight pa- 
rameters (upper and lower limits, central value and a of the 
normal function, for both T and 7). In general, the lower 
limits of f3 and 7 could be assumed to be 0.3 (strongly elon- 
gated, prolate galaxies) and 0.05 (very thin disk galaxies), 
respectively. Actually, since we are dealing with Es and cDs 
galaxies, we are allowed to explore ip down to 7 ~0.3. 

The left and right panels of Figure [3] illustrate, on 
the plane (/3,7), the result of the previously outlined bi- 
variate Lucy's rectification procedure for the WINGS sam- 
ples of Es (top panels) and BCG galaxies (bottom panels). 
For each galaxy sample, we first tried to perform the rec- 
tification without any assumption about the initial guess 
(ip° =const.). The distributions obtained in this way for the 
two galaxy samples are shown in the left panels of the figure. 
Then, we use such blind runs to refine the results by provid- 
ing more realistic initial guesses ip° , expressed as confined, 
bi-variate normal functions of 7 and T. The results of this 
refined de-projection procedure are illustrated in the right 
panels of Figure [3] For each sample, the ideal number of al- 
gorithm iterations has been choosen relying upon numerical 
simulations recording the cycle number for which the provi- 
sional and the input ip distributions display the best mutual 
agreement. Such ideal number turns out to depend on the 
sample size, being also usually greater than the correspond- 
ing values found for the one dimensional application of the 
Lucy algorithm (see Noerdlinger 1979 for perfectly oblate 
and prolate cases). 

We are obviously aware that, since (j){q) is a function 
of one variable, it cannot uniquely determine the function 
ip(P, 7). This means that the previously devised Lucy al- 
gorithm can, in principle, produce totally wrong ip distribu- 
tions which perfectly reproduce the observed distributions of 
</>(?). This actually happens for certain ip distributions (i.e. 
disk galaxies) when the first guess ip° is lacking or far from 
the true ip. In our case, extensive numerical simulations have 
shown that, even in default of the first guess (ip° =const.), 
the iterative Lucy machine is able to recover reasonably well 



the parent ip distribution. However, in order to obtain an ob- 
jective, independent check of this conclusion, we have used 
a procedure similar to that adopted by RLP93. In practice, 
assuming a gaussian (bi-variate) functional form of ip(j3,~j), 
we have scanned the (/?, 7) space with the gaussian peak, 
choosing in each point the values of the two standard de- 
viations which, according to the KS statistic, produces the 
best agreement with the observed axial ratio distribution. 
In this way, we have produced a grid of KS probabilities 
that the <p(q) resulting from the random projection of the 
local gaussian ip is drawn from the same parent population 
of the observed axis ratio distribution. The middle panels 
of Figure [3] illustrate the results of such procedure for both 
Es (top panel) and BCG galaxies (bottom panel). The im- 
pressive overlap between the regions where the distributions 
ip(f3, 7) peak according to the Lucy's algorithm (left and 
right panels) and the regions where is highest the probabil- 
ity from KS statistics (middle panels), clearly indicates that, 
in our specific case, the bi-variate version of the Lucy's al- 
gorithm works fairly well, even when it is run in blind mode 
(ip° =const.). We also note that the distribution of KS prob- 
abilities for the normal ellipticals in our sample (mid-top 
panel in Figure [3j is quite similar to the corresponding dis- 
tribution obtained by RLP93 (see Figure 5 therein). 

4.2 Results 

Figure [3] provides us with a robust indication about the in- 
trinsic shapes distributions of both Es and BCGs. Just as 
expected from the observed axial ratio distributions in Fig- 
ure they turn out to be quite apart from each other. In 
particular, there is a clear tendency of the BCGs towards 
strong prolate configurations peaked at f3 ~ 7 ~0.67, while 
the Es tend to share almost uniformely the whole range of 
triaxiality, with a slight tendency toward prolateness. 

Figure U is similar to the previous one, but it refers 
to the comparison between cD and BCG_E galaxies in our 
sample. This figure shows that the strongly prolate shape we 
found for the whole population of BCGs (see lower panels of 
Figure [3]) is entirely due to the sub-sample of cD galaxies. 
Moreover, in spite of the paucity of the BCG_E sample (re- 
sponsible for the wide distribution of the KS probabilities 
in the mid-bottom panel of Figure 3]), we note the similarity 
between the ip distributions obtained for Es and BCG_Es in 
the top-right and bottom-right panels of Figures [3] and 3J 
respectively. 

The results illustrated in Figures [3] and [4] suggest that 
cD galaxies actually constitute a foreign body inside the 
global population of cluster early-type galaxies, at least as 
far as the axial ratio distribution is concerned. To give more 
robust statistical support to this suggestion, we have ob- 
tained the axial ratio distributions of the second- and third- 
luminosity-ranked galaxies in our cluster sample and we 
have compared them with the corresponding distributions 
of Es, BCGs and cDs. The black dots in Figure [5] report, for 
the cD galaxies and for the first three ranked ellipticals in 
the WINGS clusters, the KS probabilities that the observed 
(p(q) of each galaxy sample and that of the normal Es are 
drawn from the same parent population. 

From this figure it is evident that the <p(q) distributions 
of the 2 nd - and 3 t,l -ranked galaxy samples (of which cDs just 
represent a very small fraction) are undistinguishable from 
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Figure 3. left and right panels: comparison between the probability distributions ip(/3, 7) obtained with the Lucy algorithm for the 
sample of normal ellipticals (top panels) and BCG galaxies (bottom panels) in the WINGS survey. The left panels in the figure are 
obtained in blind mode (initial guess -0° =const.), while in the right panels we let ourself be guided by the results of the blind tests 
to guess ifi. middle panels: distributions of the KS probabilities obtained, for normal ellipticals (top panel) and BCG galaxies (bottom 
panel), scanning the 7) space with a bivariate, gaussian ip and comparing in each point the expected <j>(q) with the observed axis ratio 
distribution (see text). The grey tones and the isophotes in each panel illustrate the density levels (see the uppermost grey-tone scale). 
For the left and right panels the marginal distributions of ifi onto the and 7 axes are shown in the upper and right part of each panel, 
respectively. Finally, in each panel the loci of T=0. 25, 0.5, 0.75 are reported with short-dashed lines. 



that of Es, while for the whole sample of BCGs one can ex- 
clude the hypothesis of common parent population between 
these samples and the global elliptical sample. Nevertheless, 
in analogy with the tj) distributions (see Figures [3] and [4]) , if 
we consider just the BCG_E sample (red, empty dot in the 
figure), again its <j)(l) distribution turns out to be undis- 
tinguishable from that of normal ellipticals. This confirms 
our previous claim about the peculiar intrinsic shape of cDs 
inside the family of early-type cluster galaxies. 



5 DISCUSSION 

5.1 Cluster X-ray luminosity and BCGs shapes 

The last mentioned result could also be enough to explain 
the discrepancy between our finding and the result from 
RLP93. In fact, the fraction of clusters containing at least 
one cD galaxy in our sample is ~75% (56 clusters). Instead, 
using the morphological information provided by the Nasa 
Extragalactic Database (NED), we were able to associate cD 
types to just ~26% of the BCGs in the RLP93 sample. Even 
though most of the BCGs in the RLP93 sample are actually 
lacking NED morphological information, there is an indica- 
tion that the BCG_Es (less flattened) dominate the RLP93 
sample, while the cDs (more flattened) dominate our sam- 



ple of BCGs. This circumstance can easily justify the above 
discrepancy. 

Trying to explain the remarkable difference in the 
cD frequency between the two cluster samples, we sug- 
gest it could be due to the different cluster selection cri- 
teria. In fact, WINGS clusters have been selected, in the 
redshift range 0.04-0.07, to be powerful X-ray emitters 
(LogLx >43.3 erg s -1 ). Conversely, the RLP93 sample, even 
spanning a similar range of redshift, does not obey this cri- 
terion. Therefore, the different percentages of cDs could be 
easily explained if cD galaxies preferentially resided in X-ray 
powerful clusters. Since we were not able to find in the liter- 
ature any clear indication about such a possibility, we tested 
it in our cluster sample. From Figure [6] the tendence of cDs 
to preferably reside in X-ray powerful clusters is confirmed 
with high significance (>99%) if just BCG_cDs are consid- 
ered (right panel), while the significance is slightly lower if 
we also consider nBCG_cDs (left panel). Such small differ- 
ence of significance might suggest that the above tendency is 
stronger for BCG_cDs than for nBCG_cDs. However, the KS 
test applied to the X-ray luminosity distributions of clusters 
hosting BCG_cDs and nBCG_cDs turns out to be inconclu- 
sive (Prob/cs ~50%). 

Since cDs are known to be specially luminous and size- 
able even among the BCGs, the previous result about their 
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Figure 4. Comparison between the probability distributions ip(l3, 7) for the sample of BCG_Es (bottom panels) and cD galaxies (top 
panels) in the WINGS survey. See the caption of the previous figure for details. 



preference to reside in X-ray powerful clusters suggests that 
a correlation should exist between the BCGs absolute mag- 
nitude (or size) and the cluster X-ray luminosity (or mass). 
Actually, similar c orrelations have b e en already reported 
in the literature (iBurke et all l200fi [Brough et aiTl2002l; 
Katayama et all 120031; iLin and Mohrl |2004| ; iBrough et all 
20051 ; IWhilev et alj|2008h . Figure jWshow the above corre- 
lations for the WINGS cluster sample. In this figure the 
BCGs luminosities are expressed in K-band absolute mag- 
nitudes from 2MASS (see Table [TJ and the virial masses 
of the clusters are computed from the velocity dispersion 
of galaxies inside them (both these quantities are reported 
in Table (TJ. The radii of the BCGs in the right panel of 
the figure are equivalent radii (Vox b) computed from the 
threshold area repo rted in the photometric WINGS catalogs 
l|Varela et al.|[2009h . Even if not fully homogeneous in the 
surface brightness level, this kind of size measurement turns 
out to be useful to sample the outer shapes of the BCGs. 
In Section [5.21 these shapes will be compared with those of 
dark-matter halos from ACDM A-body simulations, which 
are hardly recoverable in their inner part, due to the insuffi- 
cient resolution. Besides the Pearson correlation coefficients, 
Figure [7] reports the probabilities of the null hypothesis (no 
correlation). Moreover, the right panel of the figure also re- 
port the linear fit (r.m.s~0.13): 



Log{R k ^, G ) = 0.37(±0.056)xLo 5 (Mass c 



-3.67(±0.81)(3) 



which we use in the analysis of the ACDM A-body simula- 
tions (see next section). 



5.2 Dark-matter halos A-body simulations 

It is interesting to compare the distribution of BCG axial 
ratios in the WINGS survey with that of the central part 
(at BCGs scale) of cluster-sized dark-matter halos obtained 
with ACDM A-body simulations. In the ACDM scenario, 
small dark matter haloes form first and grow subsequently to 
larger structures via accretion and merging processes. Such 
processes are generally anisotropic, so that dark matter ha- 
los are expected to be non-spherical. A number of papers 
have been devoted to the shape analysis of dark-matter ha- 
los (see references in Section [TJ. Recently some studies, with 
both simulations and comparison with observations, have 
shown the importance of the total mass in determining the 
final halo shape. Indeed, it has been noted that the flatten- 
ing increases with halo mass, becoming more pronounced 
at Myjr > IQ 14 -^ (see, among oth ers, IWang et all 120081 : 
iMaccio et al]l2008l : lFlores et al.ll2007l ). Moreover, by using a 
set of hydrodynam i cal sim ulati ons on the cluster mas s scale , 
iKazantzidis et al.1 Hooi) and iGottlober and Yepesl (|2007n 
find that gas cooling makes the central parts of the halos 
more round. Thus, one should expect that dark matter ha- 
los of galaxy clusters are even more elongated that th e em- 
bedd ed central BCGs. Finally, both dynamical models l|Westl 
1 1994) and A-body simulations l|Dubinskilll998t ) suggest that 
halos are preferentially aligned with primordial filaments, as 
well as with the elongation of the host galaxy cluster. 

A truly realistic comparison between our observations 
and numerical simulations would require an accurate mod- 
elling of hydrodynamics, radiative cooling, star formation 
and energy feedback from SNs and AGNs, high enough mass 
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Figure 5. KS probabilities relative to the comparisons of the axial ratio distribution of Es with those of cD galaxies and of the first 
three luminosity ranked ellipticals in the WINGS cluster. The red dot refer to the BCG_E sample (non cDs). 



virial mass above 1O 14 M0, appropriate to compare to the 
WINGS galaxy clusters (see lower mass limit in right panel 
of Figure [7J . 

Given the discrete structure (particles) of the matter 
distribution in the iV-body simulations, we have used an 
iterative procedure to obtain the shapes of the halo isoden- 
sity shells. As a first guess, we calculated the main axes 
of the dark-matter particle distribution inside a sphere of 
radius Rbcg, where Rbcg is randomly assigned from equa- 
tion (3), assuming a gaussian scatter (see also right panel 
of Figure [7). Then, keeping the volume of the region fixed, 
we define a new ellipsoidal region using these axes and ob- 
tain new values for the main axes. We iterate this procedure 
until the main axes converge to a stable value, which we 
assigned to the dark matter halo as representative of its 
"simulated" BCG. The typical dark matter overdensity at 
the radii considered in this procedure is of order 10 5 times 
the background density. In the upper-left panel of Figure [5] 
we plot the distribution of the intrinsic axial ratios for the 
dark matter halo regions thus identified, 
data publicly available at http://www.mpa-garching.mpg.de/Virgo/data_d&5Moaai:dipaili son of this plot with Figures [3] and |4j led 



and force resolution, and full control over the numerical ef- 
fects. Unfortunately, our incomplete knowledge of the details 
of galaxy formation, and the insufficient computational ca- 
pabilities, make this impossible as of today. We must thus 
rely on the more robust - although simpler - modelling of the 
dark-matter distribution. In order to make such a compari- 
son still meaningful, we chose to compare the intrinsic shape 
of BCGs with that of the innermost region of cluster-sized 
dark matter halos extracted from a cosmological simulation. 

We compare the intrinsic shape of our BCGs with that 
of the innermost region of cluster-si zed dark matter haloes 
extracted from GIF2 simulation fsee lGao et al.ll2004l . , here- 
after GIF2)0, a cosmological iV-body simulation of the con- 
cordance ACDM model, performed with 400 3 dark-matter 
particles in a box of 110 Mpc/h on a side. Halos were iden- 
tified at redshift z = using a spherical overdensity cri- 
terion, cut at t he ACDM virial overdensity: A V i r = 324 
l|Eke et al.ll 19961 ). Among these we selected 3510 halos with 
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Figure 6. Left panel: comparison between the X-ray luminosity distributions of WINGS clusters hosting (thick grey line; green in the 
electronic version) and not-hosting (thin black line) at least one cD galaxy (even if not BGC). Right panel: comparison between the X-ray 
luminosity distributions of WINGS clusters with (thick grey line; green in the electronic version) or without (thin black line) BCG_cDs 
(even if not BGC). In each panel the probability that the distributions under analysis are drawn from the same parent population 
according to the KS statistics is reported. 




Figure 7. Left panel: correlation between absolute K-band magnitude (from 2MASS) of the BCGs in the WINGS survey and the cluster 
X-ray luminosity. Right panel: correlation between the virial cluster mass and the equivalent radius (in kpcs) of the BCG, derived from 
the threshold area. The error bar in the left-upper part of the plot shows the r.m.s. of the correlation (see eq. [3]and text for details) 



us to the conclusion that dark matter halos are even more 
elongated and strongly prolate than our BCGs/cDs galaxies. 
However, it is worth recalling that the axial ratios in Fig- 
ures 13 and |3] refer to the effective isophotes, while the size 
distribution we use to extract the intrinsic shapes of dark- 
matter halos in the GIF2 simulation has been calibrated 
(right panel of Figure[7]and eq.[3} on the threshold isophotes 
of the BCG sample, which, as seen in panel 15 of Figure [1] 
tend to be natter than the effective isophotes. To obtain 



a more consistent comparison, in the remaining panels of 
Figure [8] we illustrate the results of the previously outlined 
rectification procedure (with initial guess =const.) on our 
samples of BCGs (upper-right), BCG_Es (bottom-left) and 
cDs (bottom-right), when the axial ratios of the threshold 
isophotes ( q s ) are used instead of the effective axial ratios ( 
q G ). Even if in this case the BCGs tend to be more elongated 
than in the case of the effective isophote, still they turn out 
to be rounder when compared to the dark matter halos, giv- 
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Figure 8. Distribution of the intrinsic axial ratios for the sample of dark-matter halos extracted from the GIF2 TV-body simulation 
(upper-left panel), compared with the corresponding distributions (/3 s ,Tg) obtained deprojecting the distributions of the threshold 
isophotes q s for our samples of BCGs (upper-right), BCG_Es (bottom-left) and cDs (bottom-right). 



ing support to the results of iKazantzidis et al.l |2004l ) and 
iGottlober and Yepesl (|2007h about gas cooling effects on the 
central parts of simulated galaxy clusters. 



5.3 The shapes of BCGs and their evolution 

In our sample of BCGs, the distributions of many quanti- 
ties indicate that BCG_cDs are systematically different from 
BCG_Es (stronger flattening, higher mass and larger size, 
preference to reside in massive and X-ray powerful clusters; 
see Figures [5] |§] and . However, the only statistically ro- 
bust differences concern the shapes (Pks <0.001) and the 
X-ray luminosity of the host clusters (Pks <0.01). 

The tendency of cDs to reside in more X-ray power- 



ful clusters and their strong prolateness when compared to 
BCG_Es, together with the (even more) elongated and pro- 
late shapes of the inner parts of cluster-sized, simulated dark 
matter halos (see Figure [8]), could give interesting hints for 
understanding the formation histories of both the BCGs and 
the host clusters themselves. 

A detailed discussion about what the above mentioned 
results, together with those emerging from recent dynamical 
and stellar population studies of BCGs, imply about their 
formation and evolution, is beyond the scope of this paper. 
We actually plan going deeper into this topic in a forth- 
coming paper, using the whole information on BCGs from 
the WINGS database. We just provide here a few specula- 
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tions, also based on the results of some recent analyses in 
the literature. 



Stott et all l|200SF) . IWhilev et all (|200Sl ) and 



ICollins et al.l ( 200S ) have demonstrated that the stellar 



population in BCGs has been mostly in place (90% of the 
total mass) since z~2 and that, contrary to the predictions 
of cosmological si mulations and hie r archical-based galaxy 
form ation models (|Bower et al, I l2006l ; iDe Lucia and Blaizod 
2007), the BCGs are almost fully assembled a few billion 
years after the Big Bang. On the other hand, converging 
indications exist that the most luminous (and massive) 
BCGs may often display signatures of most recent star 
formation (blue-cores), likely related to the presence in the 
host clusters of cooling flows and X-ray lum inosity excess 
with respect to the averag e Lx — Tx relation (jBildfell et al.l 
120081 : lLoubser et ai1l2009h. Th e mos t luminous BCGs are 
also claimed by Bildfell et al.l (|2008l ) to be closer to the 
peaks of the cluster X-ray emission with resp ect to the 
less luminous BCGs. Finally, ICoziol et al.l (|2009l ) show that 
most BCGs are not at rest in the potential well of their 
clusters, suggesting a merging-group scenario, where BCGs 
formed first in smaller subsystems and clusters formed more 
recently from the mergers of many such groups. They also 
claim that the relative peculiar velocity inside the cluster 
(i.e. the difference between specific and average velocity, 
normalized to the velocity dispersion) is lower for BCG_cDs 
than for BCG_Es. 



All these claims, together with the findings reported in 
the present paper and the fact that the BCG_cDs are (on av- 
erage) more luminous and massive than the BCG_Es, point 
toward a scenario in which BCG_cDs and BCG_Es have ex- 
perienced quite different formation and evolution processes, 
although both have been assembled very early (old stel- 
lar populations). In particular, the BCG_cDs could have 
been formed close to the centers of quite sizeable and pro- 
late ijGottlober and Yepesl l2007h dark-matter halos, domi- 
nant with respect to the surroundig ones and progenitors 
of the most massive, present day halos (small displacement 
of the BCGs from the cluster centers and small relative pe- 
culiar velocity). In this context, considerable and prolonged 
accretion of material through cooling fl ows (enhanced X-ray 
emission) and/or dry minor mergers ((Bernard! [2009), asso- 
ciated with modest star formation (blue-cores), could have 
been at work from preferential directions (the streaming fil- 
aments), thus producing a puffing up process along them, 
leading to the formation of extended and elongated halos 
(cD morphology). On the other side, the BCG_Es could have 
formed inside group-sized halos, not much larger than the 
surroundi ng ones. Later on ; according to the scenario pro- 
posed bv lCoziol et al.l l)2009l ). the galaxy groups originating 
the BCG_Es could have merged with many other groups, to 
form cluster-sized systems. In this context, one should ex- 
pect the BCGs to show very low or absent star formation, to 
have larger displacement from the cluster centers and larger 
relative peculiar velocities compared to the BCG_cDs. More- 
over, in this case the host clusters are expected to show less 
frequently cooling flows and to be less powerful X-ray emit- 
ters. 



6 SUMMARY 

In this paper we have analysed the apparent axial ratio dis- 
tributions of the BCGs and of the normal ellipticals in our 
sample of 75 galaxy clusters from the WINGS survey. Most 
BCGs in our clusters (52) have been classified as cD galax- 
ies. The cD sample has been completed by 14 additional cDs 
(non-BCGs) we found in our clusters. We have deprojected 
the apparent axial ratios distributions of BCGs, cDs and 
normal ellipticals using a bi-variate version of the rectifica- 
tion Lucy's algorithm. Since in our case the bi-variate distri- 
bution of the intrinsic axial ratios cannot be uniquely deter- 
mined by the univariate distribution of the apparent axial 
ratios, we have used an independent Monte-Carlo technique 
to support the results of the Lucy's algorithm. Finally, we 
have used the GIF2 ACDM iV-body simulations of cluster- 
sized, dark-matter halos, to compare the intrinsic shapes of 
the inner part (BCG-sized) of halos with those of the ob- 
served BCGs. 

The main conclusions are the following: 

• The normal Es have triaxial shape, the triaxiality pa- 
rameter sharing almost evenly the (/3,7) space (with a slight 
preference for prolateness) and the 7 parameter peaking 
around 0.7. In this case, our result is fairly in agreement 
with that provided by RLP93; 

• the normal BCGs (non-cDs) and the second and third 
luminosity ranked Es in the clusters do not differ from the 
global population of Es as far as the distribution of the in- 
trinsic shapes is concerned; 

• the cDs have triaxial shape too. However, in this case 
the tendency towards prolateness is very strong and the pre- 
ferred values of /3 and 7 are significantly lower (~0.6) than 
the peak values found for Es; 

• since more than 2/3 of the BCGs in our sample are 
cD galaxies, the results of the shape analysis for our global 
sample of BCGs are similar to those obtained for cDs (pro- 
lateness, significant difference from the Es); 

• this result turns out to be strongly at variance with the 
conclusions given by RLP93, who found that BCGs and Es 
have similar apparent and intrinsic axial ratio distributions. 
This discrepancy is harldy attributable just to systematic 
differences in the axial ratios measurements; 

• since we find that, among BCGs, cDs are quite flat- 
ter that non-cDs, and since the WINGS clusters have been 
selected to be X-ray powerful, while the RLP93 clusters 
haven't, we suggest the above discrepancy to be caused by a 
preference of cDs to reside in X-ray emitting clusters. Actu- 
ally, this hypothesis turns out to be supported by the com- 
parison between the X-ray luminosity distributions of the 
cD and non-cD clusters in our sample; 

• the prolateness of the BCGs (in particular of the cDs) 
could reflect the shape of the associated dark-matter halos, 
according to the GIF2 iV-body data. 



ACKNOWLEDGMENTS 

The simulations in this paper were carried out by the 
Virgo Supercomputing Consortium using computers based 
at Computing Centre of the Max-Planck Society in Garching 
and at insitute for Computational Cosmology. The data are 
publicly available at www.mpa-garching.mpg.de/NumCos. 



Shapes of BCGs and normal Es 13 



We thank Paola Mazzei and Marisa Girardi for useful 
discussions. We are also grateful to the anonymous referee, 
whose comments allowed us to improve the paper. 



REFERENCES 

Allgood, B., Flores, R. A., Primack, J. R., Kravtsov, A. V., 

Wechsler, R. H., Faltenbacher, A., & Bullock, J. S. 2006, 

MNRAS, 367, 1781 
Bailin, J., & Steinmetz, M. 2005, ApJ, 627, 647 
Bak, J., & Statler, T. S. 2000, AJ, 120, 110 
Beers, T. C, & Geller, M. J. 1983, ApJ, 274, 491 
Bernardi M., Hyde J. B., Sheth R. K., Miller C. J., & Nichol 

R. C. 2007, AJ, 133, 1741 
Bernardi M. 2009, MNRAS, 395, 1491 
Bertin, E., & Arnouts, S. 1996, A&A, 117, 393 
Bett, P., Eke, V., Frenk, C. S., Jenkins, A., Helly, J., & 

Navarro, J. 2007, MNRAS, 376, 215 
Bildfell, C, Hoekstra, H., Babul, A., & Mahdavi, A. 2008, 

MNRAS, 389, 1637 
Bower, R. G., Benson, A. J., Malbon, R., Helly, J. G, 

Frenk, C. S., Baugh, C. M., Cole, S., & Lacey, C. G. 2006, 

MNRAS, 370, 645 
Brough, S., Collins, C. A., Burke, D. J., Mann, R. G., & 

Lynam, P. D. 2002, MNRAS, 329, L53 
Brough, S., Collins, C. A., Burke, D. J., Lynam, P. D., & 

Mann, R. G. 2005, MNRAS, 364, 1354 
Burke, D. J., Collins, C. A., & Mann, R. G. 2000, ApJLett., 

532, 105 

Carter, D. & Metcalfe, N. 1980, MNRAS, 191, 325 

Cava, A., et al. 2009, A&A, 495, 707 

Cava, A., et al. 2010, in preparation 

Cole, S., & Lacey, C. 1996, MNRAS, 281, 716 

Collins, C. A., et al. 2009, nature, 458, 603 

Conselice, C. J., 2003, ApJS, 147, 1 

Coziol, R., Andernach, H., Caretta, C. A., Alamo-Martnez, 

K. A., & Tago, E. 2009, AJ, 137, 4795 
De Lucia G., & Blaizot J. 2007, MNRAS, 375, 2 
Desroches, L. B., Quataert, E. Ma, C. P., & West, A. A. 

2007, MNRAS, 377, 402 
D'Onofrio, M., et al. 2009, in preparation 
Dubinski, J. 1998, ApJ, 502, 141 

Ebeling, H., Voges, W., Bohringer, H., Edge, A. G, Huchra, 
J. P., & Briel, U. G. 1996, MNRAS, 281, 799 

Ebeling, H., Edge, A. C, Bohringer, H., Allen, S. W., Craw- 
ford, C. S., Fabian, A. C, Voges, W., & Huchra, J. P. 1998, 
MNRAS, 301, 881 

Ebeling, H., Edge, A. G, Allen, S. W., Crawford, C. S., 
Fabian, A. C, & Huchra, J. P. 2000, MNRAS, 318, 333 

Edge, A. C, & Stewart, G. C. 1991, MNRAS, 252, 428 

Eke, V. R., Cole, S., Frenk, C. S., & Navarro, J. F. 1996, 
MNRAS, 282, 263 

Fabian, A. C. 1994, Ann.Rev.Astron. & Astrphys., 32, 277 

Fasano, G., & Vio, R. 1991, MNRAS, 249, 629 

Fasano, G., Amico, P., Bertola, F., Vio, R., & Zeilinger, W. 
W. 1993, MNRAS, 262,109 

Fasano, G., et al. 2006, A&A, 445, 805 

Fasano, G., et al. 2007, From Stars to Galaxies: Building 
the Pieces to Build Up the Universe, A. Vallenari, R. Tan- 
talo, L. Portinari, & A. Moretti, ASP Conference Series, 
374, 495 



Flores, R. A., Allgood, B., Kravtsov, A. V., Primack, J. R., 
Buote, D. A., & Bullock, J. S. 2007, MNRAS, 377, 883 

Franx, M., Illingworth, G., & de Zeeuw, T. 1991, ApJ, 383, 
112 

Gallagher, J. S. HI, & Ostriker, J. P. 1972, AJ, 77, 288 
Gao, L., White, S.D.M., Jenkins, A., Stoehr, F., & Springel, 

V. 2004, MNRAS, 355, 819 
Giacintucci, S., Venturi, T., Murgia, M., Dallacasa, D., 

Athreya, R., Bardelli, S., Mazzotta, P., & Saikia, D. J. 

2007, A&A, 476, 99 
Gottlober, S., & Yepes, G. 2007, ApJ, 664, 117 
Hashimoto, Y., Henry, J. P., & Boehringer, H. 2008, MN- 
RAS, 390, 1562 
Hoessel, J. G., Oegerle, W. R., & Schneider, D. P. 1987, 

AJ, 94, 1111 
Jones, G, & Forman, W. 1984, ApJ, 276, 38 
Katayama, H., Hayashida, K., Takahara, F., & Fujita, Y. 

2003, ApJ, 585, 687 
Kazantzidis, S., Kravtsov, A. V., Zentner, A. R., Allgood, 

B., Nagai, D., & Moore, B. 2004, ApJLett, 611, L73 
Kimm, T., & Yi, S. K. 2007, ApJ, 670, 1048 
Kormendy, J., & Djorgovski, S. 1989, Ann. Rev. Astron. & 

Astrophys., 27, 235 
Laine, S., van der Marel, R. P., Lauer, T. R., Postman, M., 

O'Dea, C. P., & Owen, F. N. 2003, AJj, 125, 478 
Lin, Y., & Mohr, J. J. 2004, ApJ, 617, 879 
Lotz, J. M., Primack, J., & Madau, P. 2004, AJ, 128, 163 
Loubser, S. I., Sanchez-Blazquez, P., Sansom, A. E., & 

Soechting, I. K. 2009, MNRAS, 398, 133 
Lucy, L. B. 1974, AJ, 79, 745 

Maccio, A. V., Dutton, A. A., & van den Bosch, F. C. 2008, 

MNRAS 391, 1940 
Mendez-Abreu, J., Aguerri, J. A. L., Corsini, E. M. & Si- 

monneau, E. 2008, A&A, 478, 353 
Merritt, D. 1985, ApJ, 289, 18 
Noerdlinger, P. D. 1979, ApJ, 234, 802 
Oegerle, W. R., & Hoessel, J. G. 1991, ApJ, 375, 15 
Oemler, A. Jr. 1976, ApJ, 209, 693 

Ostriker, J. P., & Tremaine, S. D. 1975, ApJLett., 202, 1130 
Patel, P.; Maddox, S., Pearce, F. R., Aragon-Salamanca, 

A., & Conway, E. 2006, MNRAS, 370, 851 
Paturel, G., Petit, G, Prugniel, Ph., Theureau, G., 

Rousseau, J., Brouty, M., Dubois, P., & Cambresy, L. 

1993, A&A, 412, 45 
Pignatelli, E., Fasano, G., & Cassata, P. 2006, A&A, 446, 

373 

Plionis, M., Benoist, G, Maurogordato, S., Ferrari, C, & 

Basilakos, S. 2003, ApJ, 594, 144 
Poggianti, B. M., et al. 2009, ApJLett., 697, 137 
Richstone, D. O. 1975, ApJ, 200, 535 
Richstone, D. O. 1976, ApJ, 204, 642 
Ryden, B. 1992, ApJ, 396, 445 

Ryden, B. S., Lauer, T. R., & Postman, M. 1993, ApJ, 410, 
515 

Sandage, A., & Hardy, E. 1973, ApJSupp., 183, 743 
Schneider, D. P., Gunn, J. E., & Hoessel, J. G. 1983, ApJ, 
268, 476 

Schombert, J. M. 1986, ApJSupp, 60, 603 

Schombert, J. M. 1987, ApJSupp., 64, 643 

Schombert, J. M. 1988, ApJ, 328, 475 

Schombert, J. M. 1992, Morphological and Physical Clas- 



14 G. Fasano et al. 



sification of Galaxies, G. Longo, M. Capaccioli, & G. 
Busarello, A&A, 178, 53 
Silk, J. 1976, ApJ, 208, 646 

Stott, J. P., Edge, A. C, Smith, G. P., Swinbank, A. M. & 
Ebeling, H. 2008, MNRAS, 384, 1502 
Struble, M. F. 1990, AJ, 99, 743 

Thuan, T. X., & Romanishin, W. 1981, ApJ, 248, 439 
Tremblay, B., & Merritt, D. 1996, AJ, 111, 2243 
Valentinuzzi, T., et al. 2009, |arXiv:0907.2392] 
Varela, J., et al. 2009, A&A, 497, 667 

Vanzella, E., Cristiani, S., Fontana, A., Nonino, M., 
Arnouts, S., Giallongo, E., Grazian, A., Fasano, G., 
Popesso, P., Saracco, P., & Zaggia, S. 2004, A&A, 423, 
761 

von der Linden, A., Best, P. N., Kauffmann, G., White, S. 

D. M. 2007, MNRAS, 379, 867 
Wang, Y., Yang, X., Mo, H. J., Li, C, van den Bosch, F. 

C, Fan, Z., & Chen, X. 2008, MNRAS, 385, 1511 
Warren, M. S., Quinn, P. J., Salmon, J. K., & Zurek, W. 

H. 1992, ApJ, 399, 405 
West, M. J. 1994, MNRAS, 268, 79 
Whiley, I. M., et al. 2008, MNRAS, 387, 1253 
White, S. D. M. 1976, MNRAS, 177, 717 



Table 1. Main properties of the BCG sample and of the host clusters 



Cluster 


Log(L x ) (a) 
10 44 crgs s" 1 


Log(E)( b ) 
km s _1 


z (c) 


WINGSJD^) 


Ty(«) 






JAfc 

kpc 




<$ 




(<■) 
%o 




in 1 

Is 


\°) 
<1l 


A85 


4.28 


3.022 


0.0551 


■\7i7T"\Tf~ i c Tnn/i 1 cn a k nm qi 1 c; 
WirslLrbJUU41oU.4o-Uyioll.O 


CU 


-23.89 


-26.89 


33.86 


2.5 


0.707 


0.770 


0.710 


0.640 


0.610 


0.641 


A119 


1.65 


2.936 


n n /i a a 
U.U444 


tsttmpq Tnnc;ci fi 1 o m 1 ki n n 
WllNLibJUUoolo.lz-Ulloiy.U 


cJJ 


-23.45 


-26.58 


25.67 


3.6 


0.769 


0.740 


0.690 


0.520 


n k /i k 
U.o4o 


0.767 


Aloo 


1 qo 
l.oZ 


o nno 
Z.yUo 


n n^fifi 
U.U000 


WHN LrbJUlUZ41. / Z-Z10Z00.4 


CJJ 


OQ A Q 


Ofi t^o 

-zo.oz 


OQ 7K 
Zo. / 


q n 
o.U 


U.000 


n fion 
U.OoU 


n Knn 

u.oyu 


n k en 
U.OOU 


n /i no 
U.4yo 


n K07 
U.OZ / 


A147 


0.28 


2.823 


n f\ A A 7 

U.U44 / 


tsttmp OTmnoio c\a i noi 1 qo o 
WilN LrbJUlUolZ.U4+UZlloo.Z 


rL 


-22.96 


-25.61 


26.44 


4.4 


0.769 


0.770 


0.770 


0.770 


0.789 


0.776 


A 1 CI 

A101 


n KO 

U.OZ 


Z.ool 


n okqo 
U.UOoZ 


WHN LrbJUlUoOl. lo-10Z4Zo.U 


T? 


O/i nc 
-Z4.Uo 


Ofi X.Q 

-Z0.05 


KO OQ 

OZ.Uo 


k n 
O.U 


n 77Q 

U. / to 


n onn 
U.oUU 


n 7on 
U. ( Zk) 


n 7fin 
U. f OU 


n 700 

u. /oy 


n 7nfi 
U. /UO 


A i fin 
A10U 


n 1 n 

u.iy 


o 7/i n 
z. / 4y 


n C\A QQ 
U.U4oo 


UfTMPQ Tm 1 oc;n K7 i 1 c;onoQ q 
WllNLrbJUllZOy.O/ +iozyzo.o 


CD 


00 nn 

-zz. yy 


OK Qfl 

-zo.oy 


Ofi 7Q 

ZO. 10 


Q 7 
O. i 


n 7Qn 
U. / oU 


n onn 
U.oUU 


n 7on 
U. ( Zk) 


n 7on 
U. ( ZU 


n fim 

u.oyi 


n 71 fi 
U. / 10 


A ice 
AlOo 


n Kfi 
U.00 


o 7no 
Z. /UZ 


n c\a Kn 
U.U40U 


W11N LtoJU1140 ( .Oo+UUZOOl. 1 


T? 


OQ 1 fi 

-Zo. 10 


ofi nfi 
-Z0.U0 


Ol A 7 
Z1.4 ( 


/I Q 
4.0 


n 7nn 

u. / yy 


n onn 
U.oUU 


n 7nn 
u. /UU 


n fi7n 
U.O / U 


n fion 
U.OZU 


n 7no 
U. /Uo 


Aiyo 


n 70 

u. /y 


Z.ooU 


n c\a q k 

U.U450 


\A/'TTVTr , Q TOI 0^07 fi/1 J_HS/1 1^7 
VV11N LtoJUIZOU ( .04-rUo410 I .Z 


CiJ 


O/I fi7 
-Z4.0 ( 


Ofi £.0 

-zo.oz 






n Q70 

u.o 1 y 


n ©fin 
U.00U 


n 7Qn 
U. 1 oU 


n 7fin 

U. / OU 


n 7K.0 
U. / OZ 


n kqo 
u.ooy 


A Q1 1 

Aoll 


n a i 
U.41 




n nfifii 
U.U001 


UfiMPQ Tnonnoc /ii i in/ifiQfi o 
W 11N LroJUZUyZo.41 + iy40oO.Z 


CJJ 


oq c;n 

-Zvj.oy 


Ofi CkA 

-zo.y4 


a q nn 

4o.uy 


A 7 
4. i 


n KVfi 
U.O i 


n fion 
U.OoU 


n Kon 
U.OZU 


n /ifin 
U.40U 


n A A c 
U.440 


n Kin 

u.oiy 


AO ( 


n 7i 
U. / 1 


o nQn 
Z.yoU 


n c\a ip. 
U.U4 /0 


IST'TMPQ TnO/lfinQ CkA I Qfi^/lin 1 

W 11N LroJUZ4oUv5.y4+ l 5o04iy. 1 


CJJ 


OQ c; c; 
-Zo. 00 


ofi of; 
-ZO.ZO 


a n nn 

4y.uy 


k n 

o.y 


n Q7fi 
U.o ( 


n onn 

u.oyu 


n 7nn 

u. /yu 


n 7nn 
U. /UU 


n finQ 

u.oyo 


n Qfin 

u.ooy 


A500 


0.72 


2.818 


0.0678 


iutmp Q Tn^ooxo c:i oonfion n 
Wll\ILrbJU4oooz.ol-ZzUooy.U 


cJJ 


-23.55 


-26.26 


44.62 


7 n 
( .U 


0.780 


n 7fin 
U. /oU 


0.620 


0.630 


0.662 


0.664 


A548b 


0.15 


2.928 


0.0416 


tuttmpq Tnc;/i con fio occccfi o 

WllNILTbJUo4ozy.oz-zooooo.o 


cJJ 


Ol C7 

-Zl.O ( 




4.88 


Q O 

o.z 


1.021 


n n/i n 

U.y4U 


0.790 


n 7nn 

u. /yu 


n /i 7 
U.o4/ 


0.460 


aouz 


n K7 

u.o / 


O Q K7 
Z.oO / 


n nfii n 

u.uoiy 


^s/TMr^ q Tn7coofi fii i onoi qa a 
W11N LrbJU / 0v5Zo.Ol+ZyZlo4.4 


Hi 


00 /i n 

-zz. 4y 


OK. 1 A 

-ZO. 14 


OQ Q7 

Zo.o / 


q n 
o.U 


n QQO 
U.00Z 


n onn 
U.oUU 


n 7fin 
U. f OU 


n 7fin 
U. /OU 


n 701 
U. / Zl 


n KQn 

u.osy 


A671 


0.45 


2.957 


0.0507 


WllN(^bJUoZool.OO+oUzoOo.U 


cJJ 


-24.04 


-26.72 


51.55 


5.5 


0.766 


0.790 


0.720 


0.670 


0.675 


0.698 


A 7K/1 
A / 04 


a nn 

4.uy 


q nnn 
o.UUU 


n nK/i 7 
U.U04 / 


UfTMPC TnnnoQO on nnQ7/i7 q 

w iin LrbJuyuooz.oy-uyv> / 4 / .0 


CJJ 


oq nn 

-Zo. yy 


Ofi CQ 

-ZO.Oo 


A 7 70 
4 ( . / Z 


/1 K 

4.0 


n 701 
U. / Zl 


n 7nn 
U. / UU 


n 7Qn 
U. / oU 


n fiKn 
U.OOU 


n fi k 1 
U.001 


n 77fi 
U. / / 


A VQfl 

A / oU 


QQQ 
O.OO 


O Qfifi 
Z.oOO 


U.UOvSy 


^s/TMr^QTnnicn^fic 1 ohkaq o 
W11N LrbJUyioU0.Oo-lZU04o.Z 


CJJ 


OQ fifi 

-Zo. 00 


ofi r\A 
-Z0.U4 


on 1 q 

oy. is 


/i /i 
4.4 


n QQQ 

U.000 


n oon 
U.oZU 


n onn 
U.oUU 


n 7Qn 
U. / oU 


n 7/i 
U. / 4Z 


n KQ7 

U.Oo / 


A nK7 

Ayo / 


n /i n 
U.4U 


Z.oOl 


n c\a ki 
U.U401 


UfTMPQ T1 H1 QQQ 07 OnK-KQI O 

W 11N LrbJ lUlooo.Z /-UUOOol.Z 


Hi 


oq fin 

-Zo. oy 


Ofi fiQ 

-ZO.Oo 


oq m 

zo.yi 


A Q 
4.0 


n Q77 

U.o ( / 


n o 1 n 
U.olU 


n 7nn 

u. /yu 


n 7nn 

u. /yu 


n 777 
U. / / / 


n 71 n 
U. / 1U 


A 07n 

Ay /u 


n 77 


O QQQ 
Z.OOO 


n hkgi 
u.uoy 1 


\A/'Ti\Tr~ , Q ti m 70^ 71 1 fi/ii on 

VV 11N kjoJ 1U1 1 ZO. / 1-1U4 1ZU.Z 


III 


OQ 1 fi 

-Zo. lo 


ok. /in 
-Z0.4U 






n QKfi 
U.oOO 


n 7on 
U. ( ZU 


n 7^n 
U. 1 OU 


n 7on 
U. lou 


n 7QQ 
U. / 00 


n kqo 

u.ooy 


A 1 nfin 

Aiuoy 


n a q 
U.4o 


O OQO 


n nfiKQ 
U.UOOo 


W 11N Lib J lUoy4o.44-Uo411Z.o 


Ji 


OQ fiO 

-Zo.OZ 


Ofi 1 
-Z0.1Z 


Q/1 Q7 

o4.o / 


/I K 

4.0 


n QfiQ 
U.0O0 


n 7fin 
U. / OU 


n Q7n 
U.o / U 


n Q7n 
U.o /U 


n Q7K 

U.o / 


n 7Q K 

U. / oO 


A1291 


0.22 


2.632 


0.0509 


ixnMnc ii 1 qooq oo i ccconi? n 
WllNCjbJllozzo.zz+5oooUo.U 


CJJ 


00 7n 


or 1 7 
-Zo.l / 


37.33 


A O 

4.Z 


0.745 


0.840 


0.570 


0.580 


0.513 


0.671 


A 1 fiQI ^> 

AlOola 


n Q7 
U.o / 


o Qnfi 
Z.oUO 


n c\a fii 
U.U401 


IST'TMPQ T1 OC;Q1 Q /II 1 troonQ Q 

WilN LrbJ IZO0I0.4I-IO0ZU0.0 


Hi 


OQ Q f; 
-Z0.0O 


ofi on 
-ZO.ZU 


Q1 1 K 
ol. 10 


K Q 
O.O 


n 70 q 
U. / Zo 


n 7/in 
U. / 4U 


n 7nn 
U. /UU 


n 71 n 
U. / 1U 


n fion 

u.osy 


n 7HQ 

U. / Uo 


A1644 


0.04 


3.033 


0.0467 


TTtrjMr^c ti occ7i 1 fin 1 70/IQ/1 n 
WilN LrbJ 1Z0 ( ll.oU-1 ( Z4o4.U 


CJJ 


-23.94 


-25. 14 


44.93 


2.9 


0.638 


0.630 


0.590 


0.520 


0.525 


0.557 


A 1 fifiQ 
AlOOo 


U.ol 


O Q 1 O 

Z.olZ 


n nfiQ/i 
U.U0o4 


'\X7T'\Tf~ i CT1QnQ/1fi fin I 1 Ol fi1 7 /I 

W11N LrbJ loUo4D.0U-f-iyi01 ( .4 


Hi 


OQ 07 

-Zo.Z ( 


Ofi 1 7 

-Z0.1 ( 


Ofi KQ 

ZO.Oo 


Q O 
O.Z 


n 7ki 
U. / 01 


n 7fin 
U. / OU 


n 7Qn 
U. ( oU 


n 7Qn 
U. / oU 


n 71 a 
U. / 14 


n onK 

u.oyo 


A1736 


1.21 


2.931 


0.0458 


^S/TM/" 1 C T1QOfi/1/l nn 070fi01 

WilN LrbJ lozo44.uy-z { Zozl.o 


cJJ 


-23.41 


-27.65 






0.747 


0.670 


0.630 


0.630 


0.620 


0.513 


A 1 70c; 
Al /yo 


k fi7 
0.0 / 


Z.oOU 


U.UOoo 


ISHMPQ TI £1 i OfiQCQ/1 

W UN Lrb J lo4o0Z.01+Zuo0o4.0 


CJJ 


oq on 
-Zo.oy 


Ofi fifi 

-ZO.OO 


ko nn 
OZ.UU 


/I O 
4.Z 


n 7n7 
U. / U / 


n 7fin 
U. / OU 


n 7/i n 
U. ( 4U 


n K/i n 
U.04U 


n k q 1 
U.Ool 


n finQ 
U.OUo 


A 1 QQ1 

Alool 


n n7 

u.y / 


O 7QC 

Z. / oO 


n nfiQ/i 
U.U0o4 


W11N LrbJ lo0yi0.11-(-Z / 0oo4.0 


CJJ 


O/i a n 
-Z4.4y 


07 nn 
-Z / .UU 


nc; nfi 

yo.uo 


k n 

o.y 


n con 

u.ooy 


n 7Qn 
U. / oU 


n fi Kn 
U.OOU 


n K/i n 
U.04U 


n K70 

U.O / Z 


n KQn 
U.OoU 


A 1 nQQ 


n o/i 
U.Z4 


O 700 

Z. / ZZ 


n n/i /i 7 
U.U44 / 


iutmpc Ti /i cot; c; 90 1 ifi/ioin f; 
W11N Lib J 140Z00. 00+ 104Z1U.0 


J^J 


00 fi/i 
-ZZ. 04 


OK OQ 

-ZO.Zo 


07 QQ 
Z ( .00 


/I K 

4.0 


n nQ/i 
U.yo4 


n non 

u.yzu 


n m n 

u.yiu 


n oon 
U.00U 


n o Kfi 
U.oOO 


n nKK 

u.yoo 


A 1 nm 

Aiyy i 


n fin 

u.oy 


O 777 
Z.t i t 


n nKQ/i 
U.U0o4 


W 11N Lib J 1404ol.0U-|-looooZ.o 


CJJ 


OQ QO 

-Z0.0Z 


ofi on 
-ZO.ZU 


on q c; 
0U.0O 


Q 1 

0.1 


n fiKQ 
U.OOo 


n 7nn 
U. / UU 


U.04U 


n Knn 

u.oyu 


n K71 
U.O / 1 


n fi/i 1 
U.041 


AZ1U I 


n Kfi 
U.OD 


O 770 
Z. / /Z 


n n/i i n 
U.U41U 


\s/'TTs.Tr~ i Q ti conoc no 1 01 /ific;e 1 
W11N Lib J 10oyoo.yZ+Z14u0o. 1 


T? 
Cj 


OQ A 1 

-Zo. 41 


Ofi A Q 

-Z0.4o 


O/i nK 

z4.yo 


Q 1 
O.l 


n 7Q7 

U. / / 


n o 1 n 
U.olU 


n 71 n 
U. ( 1U 


n 71 n 
U. / 1U 


n 71 k 
U. / 10 


n onK 

u.zyo 


A Ol O/I 
AZ1Z4 


n fin 

u.oy 


o nn/i 
Z.yU4 


n nfififi 
U.UOOO 


WilN Lib J I0440y.UZ+ouUOoo.y 


CJJ 


OQ 77 

-Zo. / ( 


Ofi 77 

-ZO. / / 


Qfi QQ 

0O.00 


Q O 

O.Z 


n 7Q7 

U. / / 


n 7Qn 
U. / oU 


n 7on 
U. / ZU 


n fiKn 
U.OOU 


n fion 

u.ozy 


n fi Kn 

u.ooy 


A Ol A n 
AZ14y 


n /i o 
U.4Z 


O K a Q 
Z.04o 


n nfi7n 

u.uo /y 


ti fim 00 1 1 i f;Q^fit;n q 
W11N Lib J iDUlZo.ll+OoOuOU.o 


J^j 


O/i on 

-Z4.zy 


ofi fin 
-Z0.0U 






n 7QQ 

U. 100 


n o Kn 
U.oOU 


n 7on 
U. ( oU 


n 7fin 
U. /OU 


n 71 1 
U. / 11 


n 7m 

u. /yi 


A oi fin 

Azioy 


n oq 
U.Zo 


O 7H7 
Z. /U / 


n nK7Q 
U.UO /o 


ishmpc 11 ri q^q nn i /imioo q 
W11N Lib J lolo0o.Uy-T~4yilZZ.o 


J^J 


OQ OCi 

-Zo.ZO 


OK OQ 

-ZO.Zo 






n 7/i 1 
U. / 41 


n 7on 
U. / ZU 


n 7on 
U. ( ZU 


n fifin 
U.OOU 


n fion 

u.ozy 


n 71 n 
U. / 1U 


AZZOO 


q fin 
o.OU 


o.lUO 


n nKQ 1 
U.UOol 


\S/'T'\Tr~ 1 C TI 70/107 OO 1 7QQQOR A 

W11N Lib J 1 i U4Z ( .ZZ+ ( oooZO.4 


CJJ 


oq fin 
-Zo. OU 


Ofi A A 

-Z0.44 


07 A Q 

Z / .4o 


Q K 

0.0 


n Qaa 
U.oOO 


n Q/in 
U.o4U 


n QQn 
U.00U 


n q Kn 
U.oOU 


n ono 

u.oyo 


n QQQ 

U.000 


A 0071 


n qo 
U.oZ 


o 7no 
Z. /UZ 


n nK7fi 
U.UO /0 


\A/'T'\Tr~ 1 Q TI 71 QI fi I 70m nfi 

WllNLibJl ( lolo.OD-i- ioUlUO.Z 


Hi 


OQ 7^ 

-Zo. / 


Ofi OK 

-ZO.ZO 


Kfi 1 A 

00. 14 


A Q 
4.0 


n 7nn 

u. / uy 


n 7fin 
U. / OU 


n fi7n 
U.O / U 


n fi7n 
U.O /U 


n fiQ a 
U.0o4 


n aid 
U.O /o 


A OQQO 

AZooZ 


n a fi 
U.40 


O OA Q 

z.y4o 


n nfi/i 1 
U.U041 


\S/'T'\Tr~ 1 G TOI t;i t;c; fio 1 fi/lOOl 
W UN Lib JZ10100.OZ-104ZZ1.Z 


Ji 


oq on 

-Zo.zy 


Ofi A 1 

-Z0.4 ( 


Q K QQ 

0O.00 


K 1 
0.1 


1 nKQ 
l.UOo 


n n7n 
u.y f U 


n non 

u.yzu 


n non 

u.yzu 


n QOQ 
U.0Z0 


n ki q 
U.Olo 


A OQnn 
Azoyy 


n K1 

u.oi 


O Q KO 

Z.oOZ 


n nK7Q 
U.UO /o 


isnMPC T01 C7m 70 n7cnoo n 
W11N LibJZIO ( Ul. / Z-U i OUZZ.U 


CJJ 


OO Q/1 

-ZZ.o4 


ok nfi 

-zo.yo 


1 q 1 n 

lo.iy 


Q O 
O.Z 


n 71 n 

u. / iy 


n 7/in 
U. / 4U 


n 7on 
U. (o\) 


n 7on 
U. /oU 


n O/i k 
U.o40 


n 7/i fi 
U. / 40 


A2415 


0.86 


2.843 


0.0575 


■\S7TMr^c Tooncofi 10 nc/i/iQi 1 
W11N LibJzZU0Z0.lZ-U044ol.l 


CJJ 


00 nn 


-26.14 


20.75 


5.2 


0.731 


0.720 


0.620 


0.620 


0.594 


0.592 


A2457 


0.73 


2.763 


0.0584 


WINGSJ223540.81+012905.8 


cD 


-23.53 


-26.59 


36.83 


5.2 


0.684 


0.740 


0.590 


0.580 


0.583 


0.662 


A2572a 


0.52 


2.800 


0.0390 


WINGSJ231711. 95+184204.7 


E 


-21.31 


-26.64 


3.12 


2.5 


0.901 


0.880 


0.510 


0.570 


0.577 


0.849 


A2589 


0.95 


2.912 


0.0419 


WINGSJ232357.44+164638.3 


cD 


-23.83 


-26.30 


61.08 


5.1 


0.707 


0.740 


0.570 


0.410 


0.408 


0.501 


A2593 


0.59 


2.846 


0.0417 


WINGSJ232420.08+143849.8 


cD 


-23.05 


-26.30 


21.29 


2.8 


0.650 


0.660 


0.620 


0.610 


0.613 


0.635 


A2622 


0.55 


2.843 


0.0610 


WINGSJ233501.47+272220.9 


cD 


-23.14 


-26.22 


24.05 


4.2 


0.703 


0.720 


0.650 


0.610 


0.582 


0.646 


A2626 


0.99 


2.796 


0.0548 


WINGSJ233630.49+210847.3 


cD 


-23.27 


-26.48 


25.62 


2.3 


0.666 


0.720 


0.670 


0.630 


0.618 


0.723 


A2657 


0.82 


2.581 


0.0402 


WINGSJ234457.42+091135.2 


cD 


-22.66 


-25.56 


32.60 


3.5 


0.691 


0.650 


0.630 


0.640 


0.628 


0.631 



Table 1. (continue) Main properties of the BCG sample and of the host clusters 
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(b)-(c) Dec i mal lo garithm of the velocity dispersions of cluster galaxies (b) and average redshifts of clusters (c) from lCava et all {2009). Some values have been updated, according to 
ICava et al.1 <201CD . 

(d) WINGS identifier from lVarela et al.l feOOSl) . 

(e) Morphological type (see text). 

(f) Total absolute magnitude in the V-band from lVarela et all {2009, corrected for galactic extiction) 

(g) Total absolute magnitude in the K-band from 2MASS. 

(h) Effective radius in kpc from GASPHOT. Missing values indicate that GASPHOT failed to converge. 

(i) Sersic index from GASPHOT. Missing values indicate that GASPHOT failed to converge. 

(j) Axial ratio at the effective isophote from GASPHOT. In case GASPHOT failed to converge, we derive q G from the equation: q G =0. 392+0. 572xq s (obtained from the linear 
best-fitting of the relation in panel 15 of Figure [T] 

(k)-(m) Axial ratios of the isophotes at 15 , 30 and 60 kpc of major axis [(k), (1) and (m), respectively] from GASPHOT. 
(n) Axial ratio at the threshold area from Varcla et al 1 ll2009l. SExtract or) . 
(o) Axial ratio from the LEDA Hypercat database JPaturel et"al]l2003h . 



Table 2. Main properties of the additional cD sample 
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(a) Luminosity ranking in the V-band 
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APPENDIX A: DISENTANGLING CDS FROM 
BCG_ES WITH MORPHOT 

MORPHOT is an automatic tool for galaxy morphology, 
puposely devised in the framework of the WINGS project. 
An exhaustive description of the tool will be given in a 
forthcoming paper (Fasano et al., in preparation). Here 
we just outline the logical sequence and the basic proce- 
dures of MORPHOT. It extends the classical CAS (Concen- 
tratio n /Asymmetry / clumpinesS) parameter set (jConselicel 
2003), by using 20 image-based morphological diagnostics. 
Fourteen of them have never been used, while the remaining 
six are actually already present in the literature, although 
in slightly different forms: the Sersic index, the Concentra- 
tion index dConselicel [20031 ). the Gini and M20 coefficients 
l)Lotz ct al. 200 j), th e Asymmetry and Clumpiness param- 
eters l|Conselicdl20"03h . 

Besides depending on the visual morphology, all diag- 
nostic are also empirically found to depend on the relative 
size (area enclosing 80% of the total light divided by FWHM 
area) and flattening of galaxies, as well as on the image S/N 
ratio. These 'a priori' dependencies have been removed using 
a sizeable set of simulated galaxies and a control sample of 
~1500 WINGS+SDSS galaxies, visually classified by two of 
us (GF and AD). These visual classifications have been used 
to calibrate and combine the 20 diagnostics, thus obtaining 
a final, single estimator of the morphological type, equipped 
with the proper confidence interval. This has been achieved 
by averaging the results of two different, totally independent 
approaches, based on the Neural Network (NN) and Maxi- 
mum Likelihood (ML) techniques. The first approach (NN) 
is based on a particular kind of feed-forward Neural Net- 
work, called Multila yer Perceptron Artificial NN (MLP, see 
IVanzella et~ai]|2004) . In this case the control galaxy sample 
provides the training set of the NN machine, with the mor- 
phological diagnostics used as inputs and the visual classifi- 
cations as targets. In the ML-based approach the probability 
density distributions of all diagnostics for all morphological 
types are drawn from the control galaxy sample. In this case 
the 'blind' morphological type estimate of any other galaxy 
is obtained maximizing the product of the probabilities of 
the particular set of diagnostics derived from the galaxy im- 
age. 

Defining each diagnostic and explaining its meaning is 
beyond the scope of the present paper. Here we just men- 
tion that all diagnostics are normalized to vary in the range 
(0-1) and that four out of the 14 newly devised diagnostics 
turned out to be the most effective ones in order to disen- 
tangle cDs from BCG_Es. They are: (i) a modified version of 
the Gini coefficient (GID), where the pixels are sorted ac- 
cording to the distance fro m the galaxy cen ter, rather than 
according to the flux (as in lLotz et al.ll2004h ; (ii) a Diskiness 
coefficient (DSK), measuring the correlation between az- 
imuth and pixel flux relative to the average flux value of the 
elliptical isophote passing through the pixel itself; (iii) an 
image-averaged 2D-Laplacian (LAP) in polar coordinates, 
roughly quantifying the global degree of concaveness of the 
galaxy profile; (iv) an Alignment coefficient (ALI), measur- 
ing the average concordance between the (local) maximum 
flux gradient direction and the (local) direction of the galaxy 
center. 

In Figure IAU the cumulative distributions of the four 
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Figure Al. Comparison between cDs (full line) and BCG_Es 
(dashed line; red in the electronic version) for the cumulative 
distributions of the morphological diagnostics GID, DSK, LAP 
and ALI (see text for a short description of them) . 



above mentioned diagnostics for cDs and BCG_Es are re- 
ported for comparison, together with the relative KS prob- 
abilities. The figure shows that, compared with BCG_Es, 
the cD galaxies have (on average) lower values of GID (less 
peaked profiles) and larger values of DSK (more disky), 
LAP (more concave) and ALI (local gradient better aligned 
with the galaxy center). Even if none of the four diagnostics 
illustrated in the figure, alone, turns out to be conclusive, 
their combination (through MORPHOT) provides us with 
a powerful tool for the particular task of disentangling cDs 
from BCG_Es. 



